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Abstract 

We have made an attempt to develop the quaternionic formulation of Yang 
- Mill's field equations and octonion reformulation of quantum chromo dynam- 
ics (QCD). Starting with the Lagrangian density, we have discussed the field 
equations of SU(2) and SU(3) gauge fields for both cases of global and local 
gauge symmetries. It has been shown that the three quaternion units explain 
the structure of Yang- Mill's field while the seven octonion units provide the 
consistent structure of SU (3)c gauge symmetry of quantum chromo dynamics. 
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1 Introduction 



The role of number system ( hyper complex number ) is an important factor for 
understanding the various theories of physics from macroscopic to microscopic level. 
In elementary particle physics , electromagnetism , the strong and weak nuclear 
forces are described by a combination of relativity and quantum mechanics called 
relativistic quantum field theory. The electroweak and strong interactions are de- 
scribed by the Standard Model (SM). Standard Model unifies the Glashow - Salam 
- Weinberg (GSW) electroweak theory and the quantum chromodynamics (QCD) 
theory of strong interactions. According to celebrated Hurwitz theorem [1] there 
exits four-division algebra consisting of R (real numbers), C (complex numbers), HI 
(quaternions) [21 [3] and O (octonions) [H El E]- All four algebras are alternative 
with antisymmetric associators. Real number explains will the classical Newtonian 
mechanics, complex number play an important role for the explanation beyond the 
framework of quantum theory and relativity. Quaternions are having relations with 
Pauli matrices explain non abelian gauge theory. Quaternions were very first ex- 
ample of hyper complex numbers having the significant impacts on mathematics 
& physics. Because of their beautiful and unique properties quaternions attracted 
many to study the laws of nature over the field of these numbers. Yet another com- 
plex system i.e; Octonion may play an important role [6j [71 El E] in understanding 
the physics beyond strong interaction between color degree of freedom of quarks 
and their interaction. Quaternions naturally unify [TU] electromagnetism and weak 
force, producing the electroweak SU(2) x U(l) sector of standard model. Octonions 
are used for unification programme for strong interaction with successful gauge the- 
ory of fundamental interaction i.e; octonions naturally unify [TT] electromagnetism 
and weak force producing SU(3) C X SU(2) W x U (l)y. In this paper, we have made 
an attempt to develop the quaternionic formulation of Yang - Mill's field equations 
and octonion reformulation of quantum chromo dynamics (QCD). Starting with the 
Lagrangian density, we have discussed the field equations of 577(2) and SU(3) gauge 
symmetries in terms of quaternions and octonions. It has been shown that the three 
quaternion units explain the structure of Yang- Mill's field while the seven octonion 
units provide the consistent structure of SU(3)c gauge symmetry of quantum chro- 
modynamics (QCD) as they have connected with the well known SU(3) Gellmann A 
matrices. In this case the gauge fields describe the potential and currents associated 
with the generalized fields of dyons particles carrying simultaneously the electric 
and magnetic charges. 
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2 Quaternionic Lagrangian Formulation 



Let us consider that we have two spin 1/2 fields, ip a and ip b - The Lagrangian without 
any interaction is thus defined |12j as 

L = [i^ a Yd^ a - m^M + [ iip^d^ipb - m^ b ^ b } (1) 

where m is the mass of particle, ip a and ip b are respectively used for the adjoint 
representations of ip a and ip b and the 7 matrices are defind as 

Here CTj are the well known 2x2 Pauli spin matrices. Lagrangian density (Q]) is 
thus the sum of two Lagrangians for particles a and b. We can write above equation 
more compactly by combining ip a and ip b into two component column vector; 

* = ( ) (3) 



and accordingly, there is the adjoint spinor 



i> = {^Pa i> h ) (4) 

where the spinor field ip, is the described [1 3j as a quaternion 



V>=^o + e fc V* (Vfc = 1,2,3) (5) 
followed by a multiplication rule 

ejCk = -5 jk + e jk i ei (Vj, k,l = 1,2,3). (6) 

Here 8jk and ej^i are respectively denoted as Kronecker delta symbol and three 
index Levi - Civita symbols with their usual definitions. The quaternion conjugates 
of quaternion basis elements as 

4 = e k . (7) 



3 



Accordingly the adjoint spinor if) =^70 (-0' denotes the Hermitian conjugate spinor) 
is described as 

i> = ip - e k ip k (8) 
whereas a spinor ([3]) is described as a quaternion 

V> = ipo + eii/)i+ e 2 tp 2 + e 3 ip 3 (9) 
which can be decomposed as 

ip = (ip + e^i) + e 2 {ip 2 ~ ei^s) = ipa + e 2 ip b - (10) 

This is the simpletic representation of quaternions in terms of complex number 
representations. In equation (|10|) . we have written ip a = (ip H-e^i) and ipy = (ip 2 — 
eiip 3 ) described in terms of the field of real number representations. Accordingly, 
we may write 

ip = ip - ei^i - e 2 ip 2 - e 3 ip 3 = (ip - e^i) - e 2 (4> 2 - e^) = ip\ - e 2 ipl (11) 
So, we may write the quaternionic form of the Lagrangian in terms of ^as 



(12) 



where m 



mi 



is the mass matrix with m\ is the mass of the field "01 



o m 2 

whereas m 2 is that of the field ip 2 . 



3 Quaternionic Dirac Equation 



Substituting the values of ip and ip from equations and (J3J) in equation (I12p . we 

get 
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L=i(*. *) ) < i3) 

which is reduced to equation ([1]). Defining the Eular Lagrangian equation as 



Hw)) " «0 (14) 



and taking the variation with respect to ^ a and ipb, we get 



i-f(d^ a )-mi/) a = (15) 

and 



i-yH{dM - mi(> b = 0. (16) 

Equations (I15J) and (TIB"]) are respectively recalled as the Dirac equations [13] for the 
spinors ^and ip b . Similarly if we take the variations with respect to ip a and ^ we 
get 



i{9M-f +™il>a =0 (17) 

and 



UPM-f +miJ b =0 (18) 

which are respectively recalled as the Dirac equations for the adjoint spinors ip a and 
ipb. In equations (I15H18P the 7 matrices are quaternion valued [IB] i.e. 
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1 
-1 



% e~ 



(Vj = 1,2,3; i 
These 7 matrices satisfy the following relations 



•1). 



(19) 
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(20) 



where 

9lMU = (-!,+!, +!,+!) (V/i, u = 0,1,2,3). (21) 

Let us write the Dirac equation in terms of a quaternion valued spinor ip. Now 
multiplying equation (ITtjj) by quaternion basis element e 2 , adding the resultant to 
equation (fT5|) and using equation ffTOj) , we get the Dirac equation as 

i'fid^if)) -rmp = 0. (22) 
Similarly, we may write the quaternion conjugate Dirac equation as 

i(d^)^ + m^ = 0. (23) 
Dirac equations (I22H23P provide the four current as 

f = tfrf^ (24) 
which satisfies the continuity equation = 0. 

4 Quaternionic SU(2) Global gauge symmetry 

In global gauge symmetry , the unitary transformations are independent of space 
and time. Accordingly, under SU(2) global gauge symmetry, the quaternion spinor 
ip transforms as 

ijj i y^= Uip (25) 

where U is 2 x 2 unitary matrix and satisfies 
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tfu = uW = uu- 1 = u- x u = 1. 



(26) 



On the other hand, the quaternion conjugate spinor transforms as 

i)\ — >f = i)XJ- 1 (27) 

and hence the combination ipip = ipijj = ipifj' = ijj'ip is an invariant quantity. We 
may thus write any unitary matrix as 

U = exp (i H^j (28) 

where H is Hermitian W = H. Thus, we may express the Hermitian 2x2 matrix 
in terms of four real numbers, ai^, a 3 , and 9 as 



H — 91 + o-jdj = 91 + iejdj (29) 

where 1 is the 2x2 unit matrix, aj are well known 2x2 Pauli-spin matrices and 
ex, e 2 , e 3 are the quaternion units which are connected with Pauli-spin matrices as 

e = 1; tj = —icTj- (30) 

Hence, we write the Hermitian matrix H as 



H = | 9 + a3 ai ~ m2 ). (31) 
ax + ia 2 9 — a 3 



Equation (128]) may now be reduced as 



U = exp (i 9) . exp (—ejdj) . (32) 

For £77(2) global gauge transformations both 9 and ~ct are independent of space time. 
Here exp (i9) describes the U(l) gauge transformation while the term exp (— e,-Oj) 
represents the non-Abelian SU(2) gauge transformations. Thus under global SU(2) 
gauge transformations, the Dirac spinor ip transforms as 
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■ip i — > ip ] = Utp = exp (— ejdj) ip. (33) 
The generators of this group e\ obey the commutation relation; 

[ej,e k ] = 2f jkl ei (34) 

which implies e^ej 7^ e^e* showing that the elements of the group are not commu- 
tating giving rise to the non abelian gauge structure. So, the partial derivative of 
spinor ^accordingly transforms as 

d^(x) 1 — > d^'(x) = exp (-ejdj) (d^). (35) 

As such the Lagrangian density is invariant under SU (2) global gauge transforma- 
tions i.e. 5L = 0. The Lagrangian density thus yields the continuity equation after 
taking the variations and the definitions of Euler Lagrange equations as 

^{d^j^} = ^{#7 A4 e*V'} = 9 A ,(j M ) fc = 0(VA:=l,2,3) (36) 
where the SU(2) gauge current is defined as 

{ft = {ifo^e^}. (37) 
which is the global current of the fermion field. 

5 Quaternionic SU(2) Local Gauge Symmetry 

For 577(2) local gauge transformation we may replace the unitary gauge transfor- 
mation as space- time depaendent. So replacing U by S in equation (125|) . we get 

^ ,_^'= Sip (38) 
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in which 



S = exp[-^2 qejQ(x)] (39) 
j 

where parameter £ = — J^M with ~ct(x) is infinitesimal quantity depending on space 
and time and q is described as the coupling constant. Consequently, the Lagrangian 
density (ITBl is no more invariant under SU (2) local gauge symmetry as the partial 
derivative picks an extra term i.e. 

d^^Sd^ +(d^S)^ =D^ (40) 

where the covariant derivative has been defined in terms of two Q— gauge fields 
i.e 



=d^ +A^ + B^ (41) 

where = —iA^Oj = A^ej = A^ ■ ~& and = —xB^Oj = Bfej = B^ ■ ~£ . Two 
gauge fields A^ and B^are respectively associated with electric and magnetic charges 
of dyons (i.e particles carrying the simultaneous existence of electric and magnetic 
charges). Thus the gauge field {A^} is coupled with the electric charge while the 
gauge field {B^} is coupled with the magnetic charge (i.e. magentic monopole). 
These two gauge fields are subjected by the following gauge transformations 

A'^SA.S-' + id^S- 1 ; B'^SB.S-' + id^S- 1 . (42) 

For the limiting case of infinitesimal transformations of ( , we may expand S by 
keeping only first order terms as 

S ^l + lt.-£(x); S- 1 l-lt.t(x); d M (S) (x)} . (43) 

So, on replacing partial derivative of global gauge symmetry to covariant derivative 
of local gauge symmetry, we may write the invariant Lagrangian density for the 
quaternion SU (2) gauge fields in the following form 

L = Hjry^Dnif)) - rmpip, (44) 
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which yields the following current densities of electric and magnetic charges of dyons 
i.e 

J li = (j/,) electric + Up) magnetic = ^ 1^ + *gV>7^- (45) 

where e is the electric charge and g is the magnetic charge. Equation (145 j) does not 
satisfy the usual continuity equation i.e. d^J^^O but satisfies the Noetherian form 
of continuity equation with covariant derivative as 

= 0. (46) 

6 Definition of Octonions 

An octonion x is expressed [15] as a set of eight real numbers 

7 

x = e xo + eiXi + e 2 x 2 + e 3 x 3 + e 4 x 4 + e 5 x 5 + e 6 x 6 + e 7 x 7 = e x + ^ caXa (47) 

A=l 

where eA_{A = 1,2, ...,7) are imaginary octonion units and e$ is the multiplicative 
unit element. Set of octets (eo, e±, e 2 , e%, e§, e^, e 7 ) are known as the octonion 
basis elements and satisfy the following multiplication rules 

e = 1; e e A = e A e = e A ] e A e B = -8 AB e + Jabc^c- (A B, C = 1, 2, , 7) (48) 

The structure constants f A sc is completely antisymmetric and takes the value 1 for 
following combinations, 

fABC = +1; V(ABC) = (123), (471), (257), (165), (624), (543), (736). (49) 

It is to be noted that the summation convention is used for repeated indices. Here 
the octonion algebra O is described over the algebra of real numbers having the 
vector space of dimension 8. As such we may write the following relations among 
octonion basis elements 
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[&a, e B ] = 2f ABC e c ; 
{e A , e B } = -2 5 AB e ; 

e A (e B e c ) ^ (e A e B )e c (50) 

where brackets [ , ] and { , } are used respectively for commutation and the anti 
commutation relations while 5 AB is the usual Kronecker delta-Dirac symbol. Octo- 
nion conjugate is defined as 

7 

x = e x - eix 1 - e 2 x 2 - e 3 x 3 - e^x 4 - e 5 x 5 - e 6 x 6 - e 7 x 7 = e x - ^ e A x A (51) 

A=l 

where we have used the conjugates of basis elements as eo = eo and ej = —e A . 
Hence an octonion can be decomposed in terms of its scalar (Sc(x)) and vector 
(Vec(x)) parts as 

1 1 _ 7 

Sc ( x ) = 2^ X + W ^ ] Vec ( x ) = q,( x ~~ x ^ = S 6aXa ' ^ 

A=l 

Conjugates of product of two octonions and its own are described as 

(xy) = yx; (x) = x. (53) 
while the scalar product of two octonions is defined as 

7 

(x,y) =\i,xy + yx) = \{xy + yx)= ^ x a y a . (54) 

The norm N(x) and inverse a; -1 (for a nonzero x) of an octonion are respectively 
defined as 

N(x) = xx = x x 
x- 1 



a=0 

X 

Mi) 



xx 1 = x 1 X — 1. 



(55) 



11 



The norm N(x) of an octonion x is zero if x = 0, and is always positive otherwise. 
It also satisfies the following property of normed algebra 



N(xy)= N(x)N(y)= N(y)N(x). (56) 

Equation (}5T)j) shows that octonions are not associative in nature and thus do not 
form the group in their usual form. Non - associativity of octonion algebra O 
is provided by the associator (x, y, z) = (xy)z — x(yz) Vx, y,z G O defined for 
any three octonions. If the associator is totally antisymmetric for exchanges of 
any three variables, i.e. (x,y,z) = —(z,y,x) = —(y,x,z) = —(x,z,y), then the 
algebra is called alternative. Hence, the octonion algebra is neither commutative 
nor associative but, is alternative. 



7 Gellmann A matrices 

In order to extend the symmetry from SU (2) to SU (3) we replace three Pauli spin 

matrices by eight Gellmann A matrices. Xj (j = 1, 2, 8) be the 3x3 traceless 

Hermitian matrices introduced by Gell-Mann. Their explicit forms are; 



Ai 



1 

1 




A, 



-i 

1 




A, 



1 
0-10 




1 



1 



Ab 



-i 



1 




1 
1 



A 7 






i 





—i 




As "7! 



1 
1 









(57) 
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which satisfy the following properties as 



(Aj) 1 " -Xj] 
TrXj =0; 
Tr (AjA fc ) =25 jk ; 

[Xj, X k ] = 2FjkiXi (V j, k,l = l, 2, 3, 4, 5, 6, 7, 8); (58) 

where Fj k i are the structure constants of SU(3) group defined as 



1 3 

Fl23 = 1; ^147 = ^257 = -^435 = ^651 = ^637 = ^' ^458 = ^678 = V ^ (^) 



8 Relation between Octonion and Gellmann Matri- 
ces 

Let us establish the relationship between octonion basis elements ca and Gellmann 
A matrices. Comparing equations f loD]) and (|58|) . we get 

Fabc = f abc (V ABC = 123) (60) 

and 

Fabc = ^ abc (WABC = 147,246,257,435,516,637). (61) 
Equation (160]) leads to 

= ~T~ (y A, B,C = 1, 2, 3) =► [e A , e B ] = [X A , X B ] (Ve c = iX c ). (62) 

[A A, *B\ l*C 

On the other hand equation (l6Tj) gives rise to 



=77^— (yA, B,C = 516, 624, 471, 435, 673, 572) 
[Xa, Xb\ 2iX c 

■ [e A , e B ) = [X A , X b ] (Vec = i^-)- ( 63 ) 
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We may now describe the correspondence between the matrix As and octonion units 
in the following manner i.e. 



{[e 4 , e 5 ] + [e 6 , e 7 ]} ^ --^= (e 4 e 5 - e 5 e 4 + e 6 e 7 - e 7 e 6 ) 



(64) 

Hence we may describe one to mapping (interrelationship) between octonion basis 
elements and Gellmann A matrices by using equations (I62H63P as, 



tA oc Xa =>■ £a = kXA (65) 

where k is proportionality constant depending on the different values of A i.e.k = 
i(VA = 1,2,3) and k = § (V ABC = 516,624,471,435,673,572). From equation 
(pljl . we also get 



e 3 oc A 8 =>• e 3 = kA 8 (66) 

where k = ^.With these relations between the octonion units and Gellmann A 
matrices, we may develop the octonion quantum chromo dynamics in consistent way. 
To do this, let us establish the following commutation relations for octonion basis 
elements and Gellmann A matrices i.e. 



[eg, e 5 ] [e 4 , e 7 ] e 4 1 _. 
[Afl, A 5 ] " [A 4 , A 7 ] " 2zAx " 2k, Al " ' 6lftl ' 

1 

-ie 2 fc 2 ; 
-ie 3 /c 3 ; 



[e4, e 6 ] 
[A 4 , A 6 ] 


[e 5 , e 7 ] 
[A S ,A 7 ] 


e 2 
2iA 2 


1 _^ 
2^ ^ 


• A 2 = 


[e 5 ,e 4 ] 
[A 5 ,A 4 ] 


[e 6 , e 7 ] 
[A 6 ,A 7 ] 


e 3 
2«A 3 


1 _^ 

2h ^ 


•A 3 = 


[e 7 ,ei] 


[e 6 , e 2 ] 


[e 3 , e 5 ] 


e 4 


1 


[A 7 ,Ax] 


[A 6 ,A 2 ] 


[A 3 , A 5 ] 


2iA 4 


~~ 2k A 


[e 4 , e 3 ] 


[e 7 , e 2 ] 


[ei,e 6 ] 


e 5 


1 


[A 4 ,A 3 ] 


[A 7) A 2 ] 


[Ai,A 6 ] 


2iA 5 


" 2^ 


[e5,ei] 


[e?, e 3 ] 


[e 2 , e 4 ] 


e 6 


1 


[AsA] 


[A 7 ,A 3 ] 


[A 2 ,A 4 ] 


2iA 6 


" 2h 


[ei,e 4 ] 


[e 2 , e 5 ] 


[e 3 , e 6 ] 


e 7 


1 


[Ai,A 4 ] 


[A 2 ,A 5 ] 


[A 3 ,A 6 ] 


2iA 7 


~ 2k~ 7 



A 4 = 


— ie 4 /c 4 


A 5 = 


-ie 5 k 5 


A 6 = 


-ie 6 k 6 


A 7 = 


—ie-jk-i 



(67) 

As such, we may get the following relationship between Gell Mann Amatrices and 
octonion units: 
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\ A = -ie A k A (VA= 1,2,3,4,5,6,7); 
g 

A 8 = -i-^=e 3 . (68) 

9 Octonionic Reformulation of QCD 

The local gauge theory of color Si7(3)group gives the theory of QCD. The QCD 
(quantum chromodynamics) is very close to Yang-Mills (non Abelian) gauge theory. 
The above mentioned SU (2) gauge symmetry describes the symmetry of the weak 
interactions. On the other hand, the theory of strong interactions, quantum chro- 
modynamics (QCD), is based on colour SU(3) ( namely SU(3) C ) group. This is a 
group which acts on the colour indices of quark favours described in the form of a 
basic triplet i.e. 









( R \ 






— > 


B 




\ ^3 J 







where indices R, B, and G are the three colour of quark flavours. Under SU (3) c symmetry, 
the spinor ip transforms as 



ijj i — =Uip = exp {i\ a a a (x)}ip (70) 

where A are Gellmann matrices , a = 1,2, 8 and the parameter a is space time 

dependent. We may develop accordingly the octonionic reformulation of quantum 
chromodynamics (QCD) on replacing the Gellmann A matrices by octonion ba- 
sis elements e A given by equations (|6"5]) and 066 p . Now calculating the value of 
X a a a (x) = J^o=i ^a a<1 ( x ) an d using the relations between GellMann A matrices 
and octonion units given by equations (I68p . we find 

8 

X a a a (x) = —ieikia 1 (x) — ie 2 k 2 a 2 (x) — ie 3 (k 3 a 3 (x) + k 8 a s (x)) 

0=1 

— ie 4 fc 4 a 4 (x) — ie 5 k 5 a 5 (x) — ie 6 k 6 a 6 (x) — ie 7 k 7 a 7 (x) . (71) 
Now taking following transformations 
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fcioti i — y/3 1 ; k 2 a 2 i — ► /3 2 ; (A; 3 a 3 + A; 8 
A)4Q;4 i — ► /3 4 ; k 5 a$ i — )■ /3 5 ; ^6^6 



(72) 



we get 



X a a a (x) = -ie^ 1 - ie 2 (5 2 - ie 3 /3 3 - ie A /3 4 - ie 5 /3 5 - ie 6 /3 6 - ie 7 /3 7 . (73) 

o=l 

It may also be written in the following generalized compact form i.e. 



J2^ a (x) = -tJ2e g /3 q (ar); 

9=1 



(74) 



o=l 



which may be written in terms of the following traceless Hermitian matrix form as 



9=1 



"3 + 51 



Q!i — Z«2 «4 — ^Q!5 



a 4 + za 5 a 6 + ia 7 — ^ 



(75) 



Now (169]) becomes 



$ t—Hf,* = Uif> = exp {e 9 /3 9 (x)} 



(76) 



So we may write the locally gauge invariant SU(3) C ,Lagrangian density in the 
following form; 



1 



L local = (i^D^-m^) ~-G%G^ 



(77) 



where 



= d^ + ee a A a Jj + ge a B a Jj 



(7f 



and 



g% =w - d v v; - q f abc vX c 

= (d,A a v - d v Al - e f ahc A\A c v ) + (d^B* - d v B$ - g f^B*) . (79) 
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Here in equations (ITS]) and (179|) . the e and g are the coupling constants due to 
the occurence of respectively the electric and magnetic charges on dyons. On the 
similar ground the two gauge fields {A^jand {B^} are present in the theory due to 
the occurence of respectively the electric and magnetic charges on dyons. As such, in 
the present theory we have two kinds of color gauge groups respectively associated 
with the two gauge fields of electric and magnetic charges on dyons. Hence the 
locally gauge covariant Lagrangian density is written as 

Uocai = (ifa^ - m^) - e e a A« - g (^) e a S« - \g%G? (80) 

which leads to the folowing expression for the gauge covariant current density of 
coloured dyons 

j; = e tyfij)) e a + g (V^V) e a . (81) 

which leads to the conservation of Noetherian current in octonion formulation of 
SU(3) C gauge theory of quantum chromodynamics (QCD) i.e. 

iV" = ( 82 ) 

where = J^ a X a . 
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